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NOMENCLATURE 


aspect ratio 

refers to closed circuit around a wing section 

correction accuracy 

total wing drag coefficient 

total wing lift coefficient 

total drag 

pressure drag 

skin friction drag 

spanwise force 

integral operator 

total lift 

portion of circuit C, i = 1, 2, 3, 4 
Reynolds number based on wing chord 
denotes control surface 

portion of control surface S, i = 1, . . 6 

denotes wing surface 

correction function 

denotes control volume 

free-stream velocity 

wing span 


local wing chord 
mean geometric chord 
section drag coefficient 
section lift coefficient 
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infinitesimal length vectors for line integrations 
net body force 

outward unit normal to surface S 
static pressure 
total pressure at infinity 
local total pressure 
free-stream dynamic pressure 
time 

Euler ian fluid velocity 

streainwise component of u, positive downstream 
spanwise component of u 
vertical component of u, positive upward 
Cartesian coordinates 

unit vectors along Cartesian directions x, y, z, respectively 
spanwise location of trailing portion of a horseshoe vortex 

Stt 

strength of vortex filament; circulation 
vorticity 

Cartesian components of C 
viscosity 

kinematic viscosity 
fluid density 
shear stress tensor 


Subscript e refers to an equivalent vortex 

Superscript w denotes value in region of viscous wing wake 
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DETERMINING THE LIFT AND DRAG DISTRIBUTIONS ON A 


THREE-DIMENSIONAL AIRFOIL FROM FLOW-FIELD 
VELOCITY SURVEYS 
Kenneth L. Orloff 
Ames Research Center 


SUMMARY 


The application of the incompressible momentum integral equation to a 
three-dimensional airfoil is reviewed. The objective is to interpret the 
resulting equations in a way that suggests a reasonable experimental techniqu 
for determining the spanwise distributions of lift and drag, Consideration 
is given to constraints that must be placed on the character of the vortex 
wake structure shed by the wing, to provide the familiar relationship between 
lift and bound vortlcity. It is shown that the induced drag distribution is 
not directly measurable, but can be obtained, via the lift distribution, 
approximately for a. deflected wake and exactly for a planar wake. A novel 
technique is presented for obtaining the spanwise lift distribution from 
velocity surveys behind the wing. Moreover, it is shown that it is only 
necessary to survey a short distance above and below the wing trailing edge. 
While the measured lift coefficient is not the true value, it is, neverthe- 
less, accurate. The necessary formalism is developed to correct these 
measured values by using an equivalent single vortex model to account for the 
unmeasured portion of the downward (or upward) momentum. Examples are 
presented for several typical loading distributions and the results of a 
numerical simulation of the suggested experiment are discussed. 


INTRODUCTION 


This paper deals with the application of the momentum integral equation 
to an airfoil that is generating both lift and drag. The objective is to 
interpret the resulting equations in a way that suggests a reasonable 
experimental technique for determining the spanwise distributions of lift 
and drag on a three-dimensional wing. 

The motivation for the analysis has been largely due to the development 
and success of the laser velocimeter for conducting wind-tunnel flow diagnos- 
tics. Conventional velocity measuring techniques are either cumbersome or 
lack the required accuracy for conducting flow surveys. With the laser 
velocimeter, however, these drawbacks are overcome, and one may assume (as 
is done in the momentum analysis to follow) that accurate three-dimensional 
velocity data are available, as may be required in the application of the 
resulting equations. 






The guiding philosophy in the analysis herein is that one should be 
dutifully aware of the true meaning of the data that are obtained experimen- 
tally. Flow characteristics that are easily measurable should be measured 
accurately; those that are not directly measurable with sufficient accuracy 
should be either replaced by equivalent measurements that are accurate or 
should be computed from an equivalent model of the flow, the parameters of 
which are based on the experimental data whose limitations are well understood 


MOMENTUM ANALYSIS 


The integral form of the momentum equation may be written as 

J (pu)dV + ^ pu(n • n)dS = ^ pf dV - ^ pn dS + ^ (n • ^)dS (1) 

V S V S *^S 

-y 

where f is the net body force per unit mass at points within the volume 
V, and T is the shear stress tensor at points on the closed surface S 
surrounding V. As is conventional, n is the outward unit normal^to the 
surface S, p is the static pressure, p the fluid density, and u the 
Eulerian fluid velocity. If the flow is assumed to be steady, and no body 
forces are present, then equation ( 1 ) simplifies to 

^ pu(u • n)dS = - ^ pn dS + ^ (n • ^)dS (2) 

s s s 

Two different control surfaces S are considered. Figure 1 depicts a 
control surface that would be used to compute the total lift and drag on the 
airfoil. The effects of the airfoil on the fluid are represented by corre- 
sponding variations in the velocities, static pressures, and shearing 
stresses on the control surface. 

To compute the sectional lift and drag, the control surface in 
figure 2 is used. This surface is basically the same as that of figure 1, 
except for the infinitesimal thickness Ay. Also, while the wing surface 
forms part of the closed surface S, the portions on S 5 and Sg in figure 2 
which cover the airfoil cross section are not part of the closed cpntrol 
surface. 

For either of these control surfaces, the momentum equation (2) may be 
written as 
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- f puu dy dz + J pCiu dx dy + J puu dy dz - T puu dx dy 

^1 ^2 ^3 

- J puu dx dz + j puu dx dz + J pu(u • n)dS 

*s. y ^ S 

5 6 w 

- - f pH dy dz - f pn dx dy - f pn dy dz “ f pn dx dy ~ J* pn dx dz 

T-i ^r% ^r% 


pn dx dz - J pn dS + 


^ (n • ^)dS + ^ (n • t)dS 


thru Sg 


We note the following four observations with respect to equation (3): 

• J pu(u • n)dS vanishes since u • n = 0 at all points on the 

S 

w 

wing surface. 

• Appendix A shows that f (n • T)dS =0 for i = 1, . . 6 

^i_ 

when the Reynolds number (cV^/v) » 1. 

• J (n • t)dS = “D ^x where D ^ is the drag due to skin friction. 

•'S sf sf 

w 

Also, we assume the lift generated by the shearing forces at the 
surface of the wing to be negligible compared to the pressure lift; 
hence, this integral has a negligible z-component. 

• - r pn dS represents the net pressure force exerted by the airfoil 

S 

. w , 

surface on the fluid. This is exactly minus the force exerted by 
the fluid on the airfoil. In component form, 


- f pn dS = ~D X and 
•'s ^ P 


- / pn dS = — Lz 
w 


where both the form drag and the induced drag are contained in D^. 

Incorporating these observations into equation (3), we obtain, after 
separating into components , 





L = I pu u dy dz - / (pu ^ + p)dx dy - I PU u dy dz 
Jc ^ ^ Jc ^ ^ ^ 




+ ^ (pu_^ + p)dx dy + / dx dz 
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pu u dx dz 

y z , 


= D _ + D_ = / (pu„2 + p)dy dz - / pu u dx dy - / (pu ^ + p)dy dz 

sf P X ^ ^ •'^3 

+ J pu_u_ dx dy + f pu^u^^ dx dz - ^ pu^u^^ dx dz 


X z 
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F = I pu u dy dz - I pu u dx dy 


X y 


-X 
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X. 

- / pu u 

M ^3 
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j pu u dx dy + I (pu,^ + p)dx dz - f (pu ^ 4- p)dx dz 

« y z Jr y Js y 
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Fy, a spanwise force, can exist for the control surface of figure 2, but 
vanishes for a wing having bilateral symmetry with a control surface as 
shown in figure 1. 

Several of the terms in equations (4) require a knowledge of the static 
pressure on the control surface. Experimentally, these pressures would be 
difficult to measure accurately; undetected, small changes over a large 
surface area can produce a considerable error in the resulting force computa- 
tion. Since the accurate measurement of the velocity is somewhat more 
convenient (especially with the laser velocimeter) , the static pressures 
in equations (4) are replaced with the total pressure and the local velocity 


P = Pt - I V + Uy^ + u/) 

This substitution is helpful in two ways: (1) a difficult measurement has 
been replaced by simpler ones; (2) the total pressure p^ is the total head 
Pq in all regions that have not suffered viscous losses. 

Using this replacement of the static pressure and nondimensionalizing 
the velocities by and x, y, and z by c, the mean geometric chord, 

we obtain 
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— = 2f uudydz-f (u^-u^ - u 2)dx dy - 2 /* u u dy d: 
J„ X z Jc, z X y X z 


•'S 


+ r (u^-u^-u 2)dx dy + zf u u dx dz - 2( u u dx dz 
Jo ^ X y •'S ^ •'^6 y 

^4 5 

— = r (u ^ - u 2 - u ^)dy dz - 2 f u u dx dy 

^co •'s. ^ y ^ ^ ^ 


f (u2_u2-u 2)d:y dz + 2 f u u dx dy + 2 J u u 

X y z Jg X z •'S ^ y 


dx dz 


^ 5 ) 


/. 


?00 


dy dz 


- 2 I u u dx dz + 

Jo y 

— = iC u u dy dz - 2/* u u dx dy - 2/* u u dy dz + 2 ( u u dx dy 

q„ Js^y •'s y ^ •'s, y •'s y ^ 

1 2 - ■J 

+f (u^-u^-u ^)dx dz - I (u 2 - u 2 _ u 2)dx dz 
J y x z Jgy x z 

^ ^5 6 

where p is the total pressure at infinity and p^ is the local total 
pressure at any point on S . 


Total Lift and Drag 


We now proceed to let all surfaces of the control volume in figure 1 
expand to infinity, with the exception of surface S3. To maintain Fy = 0 , 
S5 and Sg expand symmetrically about the centerline of the wing. On these 
remotely located surfaces, Uy ^ u^ 0 and ^ !• if we substitute into 
equations ( 5 ) the coefficient representations for lift and drag, L = q^CL(bc) 


and 





> ( 6 ) 


dy dz 


J 


where the integration is over the rear surface S^. 


In equations (6), the analytical expressions for the lift and drag 
coefficients require a knowledge of the velocity and total pressure distribu- 
tions in the wake of the airfoil only. On the other hand, it is unlikely 
that the implied technique would be implemented experimentally, since the 
region of the wake that must be measured is extensive; further, the procedure 
would be influenced by the interference of the wind-tunnel walls. Even if the 
entire flow could be mapped accurately, the same values could be obtained with 
much less difficulty by using a force balance. Of greater interest are the 
wing sectional properties, which cannot be obtained from force measurements. 


Section Lift and Drag 


We return to equations (5) and consider the control surface presented in 
figure 2. The differential dy is now Ay and the area integrations over 
Sj, S^, Sg, and become line integrations along Lj , L^, Lg and 
respectively, with width Ay. L and D are represented by L = q^c^(cAy), 
and D = q^oC^jCcAy), where the coefficients, and c^, are now the section 

values and c is the local chord. Equations (5) are now written as 


= 2 


if u u 1 dz - / 

^ Z Jt 


(u ^ - u ^ - u ^)dx - 2 1 u u dz 


f. 


X z 


dx dz - 




u u dx dz 

y z 
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+ r (u^-u^-u ^)dx + “ f f u u 
•'l ^ X y Ay y 

^ 5 ; b 

= r (u^“U^-u ^)dz - 2f u u dx - f (u ^ - u ^ - u dz 
a Jj X y z J X z X y z ' 

If uu dx dz - I uu 
Us X y J. X 


} (7) 


“1 
+ 2 


u u dx + - — 
^ X z Ay 

^4 


J. 


I (^) 


dx dz 




dz 


The equation for the spanwise force, Fy, has been dropped from the analysis 
at this point because it is found not to lead to anything of special interest 
for configurations that have bilateral symmetry. 


The assumption is now made that the flow velocity around the airfoil is 
changing smoothly and slowly enough with respect to y that the integrands 
in the remaining surfaces integrals in equations (7) can be expressed to 
first order as 


(u u )_ 
y z 6 


(u u ) . 
X y 6 


(u u ) ^ + 
y z 5 

[i 

Ay 

s 

(u u ) + 

(u U )] 

Ay 

X y 5 

[ay xyj 

S ^ 



where subscripts 5 and 6 refer to the surfaces Sg and Sg, respectively, 
and the subscript S indicates evaluation of the derivative at the surface 
spanned by the closed circuit defined by , L^, Lg, and L^^. The normal 
to S is taken in the positive y~direction. Spbstituting equations (8) into 
the surface integrations of equations (7) , 


6 
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( u u dx dz - I u u 

o y 2 Jq y 2 


dx dz 


f u u dx dz) = -2 I 

■'s y ^ / Js 

b 

I u u dx dz I = -2 I 

•'s. ^ y / Js 


(u u )dx dz 


3y y z 


(u u )dx dz 


9y X y 


Appendix B shows that the surface integrals on the righthand side of equa- 
tions (9) can be written as 

-2 1 (u u )dx dz = -2 r u u dz + j (u ^ - u ^ - u ^)dx 
Js 9y y 2 X z z X y 

+ 2| u u dz - r (u ^ - u ^ - u 
J X z Jj z X y 

/•/• /8u 3u \ /3u 3u \ 

*'ffs '‘y 


-r~ (u u )dx dz 
3y y z 


= - f (u^-u^-u ^)dz + 2 j u u dx 

•'L ^ ^ ^ •'L ^ ^ 

1 2 

+ / (u^-u^-u ^)dz - 2 | u u dx 

J X y z J X z 


£ ( 


3u 3u 
__X _ 

3x 3y 


dx dz 


-flM 


3u 3u ' 

X z 

3z 3x 


dx dz 


Equations (7) may, therefore, be reduced using equations (10) 


cc. ~ 2 fC (u ^ - u C )dx dz 

a JJ^ X y y X 

cc, - zff (u C - u 4' ''dx dz + f 

d JJs y 2 ^ y Jl, V 7 


where C = curl u is the vorticity vector. 

Equations (11) are not yet in a form to suggest a reasonable experimental 
technique for determining the aerodynamic coefficients Cjj, and c^^. The method 
implied by these equations requires detailed distributions of the velocity 
field over the entire surface S. The additional need for the vorticity 
distributions requires that the derivatives of the velocity be determined from 
the data, reducing the accuracy of coefficient results* Hence, before our 
analysis can continue, some assumptions must be made regarding the character 
of the flow, in order to further simplifty equations (11) . 
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Lifting-Line Approximation 



If the flow field of interest can be adequately described by a lifting- 
line model, at least over the region enclosed by the closed circuit C, then 
equations (11) can be simplified considerably. The assumed geometry, shown in 
figure 3, reduces the wing cross-sectional area to zero about the bound 
vortex. The surface denoted by the double integration in equations (11) is 
therefore the entire area spanned by C. This allows the use of Stoke’ s 
Theorem, which reduces the problem to a line integral around C, if equa- 
tions (11) can be put into a form that is compatible with this theorem. 


To apply Stoke’ s Theorem, we note two implicit restrictions of the 
lifting-line model: 

• If the lines of shed vorticity are allowed downward deflection, 

due to mutual Induction, but spanWise deflection is forbidden, then 
Cy is nonzero only along the bound vortex. Hence, it is only 
nonzero at x = z = 0, where we must have u^^ = 1. 

• Cx vanishes along the bound vortex. It is nonzero only along a 
trailing line of shed vorticity of strength dr(y). However, if 
no spanwise wake deflection is allowed, then Uy = 0 qlong these 
lines . 

Using the two restrictions above, the lift coefficient is 


cc- = iff (u 5 
a x^y 


u C )dx dz = 2 ff C dx dz 

y x' y 


Writing Cy = y • curl u we may apply Stoke’ s Theorem, 

^^0 “ ^Xf y * ^ *^x dz = 2® u • dil (12) 

^ " •'C 

where y is precisely the unit normal to the surface S. This is a familiar 
result for two-dimensional flows, but it is also valid for a three-dimensional 
flow with the two lifting-line restrictions above. 


If the lifting-line model is used to simplify the drag coefficient, the 
"planar wake" assumption must be made, which additionally requires 
whereas this value was not restricted for the lift coefficient. With this 
additional restriction, ^ ^ 




dx dz + 



Since Cv is nonzero only at x = z = 0 on S, then 


cc 




dz 


i 


1. 


where u^, is the value of the downwash at the bound vortex at the point 
where it is intersected by S. This Is further reduced to 


CCd = -u^(cc^) + 



Pn - Pr 


dz 


(13) 


The first term in equation (13) is precisely the induced drag. However, 
it is not possible to measure experimentally the induced downwash at the 

chordwise location of the wing center of pressure (bound vortex). Hence, from 
lifting-line theory we replace Ug with (in dimensional form) , 


u^(y) 


J b/2 3r(v') 

3y' 

y - y' 

-b/2 


dy’ 


(14) 


Writing P(y’) “ (l/2)Vooc(y’ )c^(y’ ) , the nondimensional form of equation (14) 
becomes . 

►b/2c 

^ I dy> (15) 


Sir 


J b/ zc 
h/z'^ 


y - y 


Substituting equation (15) into equation (13) , 

b/2E 


CCj = 

a 


8.J_ 


b/2c 


3/3y' (cc^) 

y - y' 


dy ’ + 


I 


Pn “ Pr 


dz 


(16) 


Hence, the induced drag is not explicitly determined, but is only implicit 
via the lift distribution under severe restrictions (i.e., no vertical or 
spanwise wake deflection). At any rate, it is first necessary to determine 
the lift distribution accurately from equation (12). 


EXPERIMENTAL CONSIDERATIONS 


If the lifting-line model is assumed to represent the real flow ade- 
quately, at least in the region of the control surface, then equation (2) 
suggests that one can determine the section-lift coefficient experimentally 
by measuring the appropriate velocity components on the closed contour C. 

If the orientation of the loop remains rectilinear to the flow, then the 
extent of circuit C is not restricted as long as the vorticity passing 
through it remains constant. As the loop expands, however, the velocity vari- 
ation about free-stream conditions becomes smaller; to sense these changes 
with sufficient accuracy to perform the integration in equation (12) becomes 
increasingly difficult. 







To measure the flow velocities along the loop C with a laser veloclai’^ 
eter, one must have "spanwise" optical access to the flow, as indicated in 
figure 4. Such a view of the selected spanwise location may, however, be 
hindered by obstacles such as engine nacelles, flap pods, etc. To measure 
around these obstacles, the loop may become so large that the accuracy of the 
measurement is compromised. Additionally, spanwise access may not be feasible 
as in a large wind tunnel (e.g. , Ames 40- by 80-Foot Wind Tunnel) . See 
figure 5. Hence, it would be convenient to extend the results of the previous 
sections to the spanwise loading when only limited optical access is available 
The next section presents such an extension of the theory. 


Lift Distribution from Wake Surveys 


Equation (12) can be written 




If the boundaries Lj , Lp, and are allowed to expand to infinity, 
u„ = 0 on L, and u„ = 1 on Lo and Lh. Then we have 

Z 1 X H 



then 


(17) 


Because an experimental traverse along Lg can be made over only a finite 
distance between and z^ (see fig. 6), we express equation (17) as 

cc^ = -2 I u^ dz - 2 I u^ dz + (cc^) (18) 

Jz^ J -00 meas 

, I ■ 

where > 

A2 

(cc^) = -21 u dz 

meas •'z^ 

is the "measured” value between the limits Zj and The first and second 

terms in equation (18) relate to the z-component of the linear momentum, 
which is present outside the measured region. This momentum must be accounted 
for if we are to have an accurate value for the section lift coefficient. 


The difference to be expected between the true value of cc^ and the 
measured value (^^ 0 } can be shown by formulating a computer code to 

meas 


Because the experiment is simulated precisely on a minicomputer , 
the numerical computation of the lift coefficient is referred to as the 
"measured" value; the quotation marks are intended to distinguish this 
from an actual wind-tunnel measurement. 
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generate a flow field and then performing the experiment numerically. 

Appendix C presents the details of this flow- field generation routine and the 
* relevant equations. Note that the generator allows for a nonplanar wake, but 

there is no spanwise vorticity, ?y, in the free wake. Figure 7 shows, in 
flow chart form, the experimental simulation that uses this velocity field as 
a basis for the "measurements." This routine is carried out on a high-speed 
minicomputer with a graphics display terminal. The figures shown in Appen- 
dix D are reproductions of the hardcopy from this terminal. 

Figure 8 compares the known ccj^ distribution with that determined by a 
wake survey over a finite distance, z-^ to The points for the "measured" 

loading have been obtained from the experimental simulation routine. The 
flow-field generator produces a superposition of two horseshoe systems, one 
shed inboard at y = 1.32 and the other near the tip at y = 3.38. The core 
radius of the outer vortex, set to 0.12, provides ccn - 0 at the wing tip, 

; y = 3.50, The survey limits have been chosen as one mean geometric chord c 

f both above and below the "trailing edge"; the survey is made at 0.9 mean chord 

behind the bound vortex. The comparison is understandably poor because a 
^ significant portion of downward momentum is contained in the unmeasured 

: region and has been neglected. Also, note that the "measured" loading in 

; figure 8 does not extend inboard from y =0.5. This has been done inten- 

tionally since it is unlikely that one could obtain data any closer in a real 
U experiment because of fuselage width. 



Correction of "Measured" Lift Distribution Using Equivalent Vortex Model 

The difference between the "actual" and "measured" loadings in figure 8 
can be accounted for by modeling the flow analytically in the unmeasured 
regions. The parameters of the model are determined by the character of the 
data obtained in the measured region. The correction terms in equation (18) 
can then be computed and the corrected lift coefficient determined. 

The assumption is made that there exists a single "equivalent" horseshoe 
vortex that contains an amount of downward momentum in the unmeasured regions 
that is very nearly equal to the momentum generated in these regions by the 
span loading of interest. If it is further assumed that this horseshoe is 
planar, then the induced velocity component u^ at any point (x, y, z) in the 
flow due to the equivalent system can be shown to be (see appendix C), 



where the subscript e refers to the equivalent vortex, = (cc^) /8ir is 
the strength of the vortex, y^ is the spanwise location of the trailing 
portion of the horseshoe, and A, B, and C are given by 
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= r 1 

|_»^ + + y ^ + z^ J 

® " f'r^=%=T'^ 

.l_/x- + yj- + z2 J 


(x^ + z^) 


^2 ^ y 2 ^ 2^ 

+ 


2 a. ,2^-l 


+ 1 


where y^. = Yq + y and y_ = Yq “ y. 

e e 

Substituting equation (19) into equation (18), 

cc = 2y [xI(A) + y_I(B) + 74 . 1 ( 0 )] + (cc^) 


where 1(F) is the definite integral operator 


1 (F) =1 F dz + I F dz 


The indefinite form of the integrals in equation (21) may be evaluated (and 
shown to be vi 


/* - 4 [»- (^) ^ (^^)] 


J ' ^ *. -II X s 

B d z = — ■ tan I 

y_ V y 


sin 6 - 


+ y ^ + z^ 
4- 


and sin ^ = 


i4c^ + V ^ + z^ 


With attention being given to the signs of the numerators and denominators 
in equation ( 22 ), equation ( 20 ) can be expressed as 



2tt 


(cCjj)^T^ + (ccj,) 


-Yn s y sy, 




for y > Y( 




y < -yo 


where 


/y^ sin _ 

T^ = tan ■‘- { — — 1 + tan ^ 


X sin 0, 


± tan" 


y sin 


± tan*"^ 


X sin c}>^Y _i/^2\ _i 

I + tan ( — I ± tan 

y- / \y^J 


m- 


7T + IT 


This solution is restricted to Z 2 = -z-^; subscript 2 in equation (24) 
denotes evaluation of sin 0 and sin c|) at z = Z 2 * '^^e distance behind the 
bound vortex at which the survey has been made is given by x, and y is the 
spanwise location where the correction is desired. 

Note, from equations (23) and (24) , that the correction is completely 

determined by specifying two values: (tc ) and y„ , the strength and loca- 

£ e 

tion, respectively, of the equivalent vortex. To obtain these values, two 
known conditions are invoked relevant to the corrected distribution 
These two conditions are: 

• The value of cCj^ at the wing tip must be zero. 

• The area beneath the cc£ distribution must equal the total lift 
coefficient of the wing, C^. 

The first condition is satisfied only if y^ is located inboard of the wing 
tip. It then follows from the second equation of equations (23) that 

at tip (25! 


T_ contains both y^ and (cc^)^ so that closure of the solution requires 

e ; ■ 

the second condition to be invoked. This is accomplished numerically by 
iterating the solution, as indicated in figure 9. 
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Results from Numerical Experiments 

Figure 10 shows the results obtained by applying this iteration scheme 

(fig. 9) to the data in figure 8. The equivalent vortex is located at 

y« = 2.79, with a strength (cc ) = 1.A2. The area beneath the corrected 

^ e X, e 

data matches = 1.08 to within 1%. 


The agreement between the corrected data and the "actual" loading in 
figure 10 is good, except in the areas noted. If, however, the extent of the 
measured region were greater than - 1,0, in this case, then the correction 
scheme should be more accurate. To verify this, we refer to figure 11, which 
shows uncorrected loadings for several survey distances from 0.5 to 

2.0 chords above and below the wing. Clearly, as survey distance 
increases, more momentum contributes to the "measured" value. Figure 12 
presents the corrected loadings for the most extreme values of z^^ in fig- 
ure 11. As expected, the loading is more faithfully reproduced at z,^ ~ 2.0 
than at 5. 


To quantify the accuracy with which this correction scheme reproduces 
the "actual" loading, the area mismatch in figure 10 is used and the following 
definition is applied; 


Correction accuracy 




area mismatch 


)■ 


100 % 


(26) 


The correction accuracy (abbreviated CA) in figure 12 indicates a mismatch 
of less than 4% for the survey of 2.0 chords, compared to more than 20% mis- 
match for a survey of 0.5 chords. 


It is of interest to consider the application of this equivalent vortex 
correction technique to several different loading distributions. The loadings 
chosen to represent typical aircraft configurations are described in table 1. 
The results obtained for wake surveys of these loadings at x = 0.9 for 
Zg, = 0.5 and ~ 2.0 are presented in figures 13, 14, and 15. Interipediate 
values of z^ and a complete accounting of the results are provided in 
appendix D. 


Figure 15 indicates that, for the elliptic distribution, a correction 
accuracy of 95. 8% can be attained for a wake survey of only 0,5 chords above 
and below the wing. This finding is in sharp contrast to the value of 78.1% 
obtained for the flap/tip combination (fig. 12) with the same survey distance. 
Moreover, the 30®/ 0° and 30°/30° loadings indicate correction accuracies that 
are intermediate between these values. To understand this trend more thor- 
oughly, we refer to figure 16, which shows the correction accuracy for the 
loadings studied as a function of survey distance z^^^ These data suggest 
that as the loading distribution deviates further from the elliptic, the 
required survey distance must increase to maintain a given percentage accuracy 
in the agreement with the "actual" distribution. Equivalently, this means 
that the single vortex equivalent model becomes less adequate as the loading 
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shifts predominantly inboard, exhibiting less resemblance to an elliptic load- 
ing. On the other hand, appendix D reveals that an acceptable reproduction of 
the loading is obtained with the single vortex correction when the correction 
accuracy is greater than about 93%. Therefore, if one were to conduct the 
wind-tunnel experiment with a laser velocimeter, a minimum survey distance 
would be required to obtain an accurate (93%) loading distribution. Figure 16 
suggests the following guidelines for any loading similar to those presented 
herein: 


Loading 

Flap /tip 
30/0 


Minimum survey distance 
(93% accuracy) 

1.4 

1.0 


30/30 


0.8 


Elliptic 


0.3 


Distribution of Induced Drag 


Equation (16) provides a means for determining the spanwise distribution 
of induced drag: 





3/Sy'(ccj) 

y - y' 



(27) 


This expression was obtained by assuming that: (1) the wake is entirely 

planar (Cy = K.z ~ y (2) the induced drag can be expressed as the product 

-U2,(ccj^) with Uz obtained from the lift distribution. Even though the flows 
in the previous section have induced downward deflection, it is instructive 
to assume that the "actual" induced drag can be represented by -U2(cc^), and 
that the computed drag can be obtained numerically from equation (2'7) using 
the cc£ data resulting from the wake surveys of the previous section. 

Figures 17, 18, and 19 present the induced drag distributions obtained 
from the loadings in figures D-9, D-14 , and D- 20 of appendix D, for the 30/0, 
30/30, and elliptic cases, respectively. As the loading shifts predominantly 
inboard (30/0 and 30/30), the high induced drag at the tip decreases, as 
expected. These figures must, however, only be used to infer regions of high 
and low induced drag and overall trends, since equation (27) contains assump- 
tions that cannot be closely satisfied. Even the "actual" drag shown is not 
precise because the wake is not truly planar, and -^z is obtained from the 
known (deflected) flowfield. Nevertheless, the three examples in these 
figures have been included for completeness. 


Conveniently, the velocity wake surveys may be conducted along the same 
path L3 as that used to determine the viscous drag. Knowledge of the three- 
dimensional velocity along this line allows one to preset the orientation of 
the total pressure probe so as to align it more precisely with the mean flow 
direction, thereby improving its accuracy. 
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SUMMARY AND LIMITATIONS 





An analysis of the momentum integral equation has shown that the local 
section-lift coefficient on a wing is determined solely by the net vorticity 
passing through a closed circuit C when certain restrictions are placed on 
the character of the inviscid wake structure. It was shown that, when downward 
deflection of the shed vorticity is present, but spanwise deflection is absent 
within the circuit C, the relationship in equation (12) is valid. 

The analysis has assumed the bound vortex line to lie along the y-axis 
with analysis does, however, proceed identically from equa- 

tion (11) when wing sweep is considered. In this case, Cx has a finite 
value along the bound vortex, but spanwise deflection of the trailing vortex 
system is again forbiddenj and symmetry provides Uy = 0 at all points on the 
bound vortex. Hence, the examples presented for the swept-wing transport-type 
loadings are still meaningful. 

A method has been developed for determining the spanwise lift distribu- 
tion from wake measurements only. The equivalent vortex correction technique 
has been shown to be adequate as long as the wake survey line extends suffi- 
ciently far above and below the wing trailing edge. The required survey 
distance has been shown to be governed by the degree to which the loading is 
concentrated inboard. 



TABLE 1.- LOADINGS USED EOR FLOW-FIELD SIMULATION 


Loading 

designation 

Number of horseshoe 
vortices used for 
flow-field generator 

^L 

Remarks 

Flap/tip 

2 

1.08 

Loading presented in figures 8, 
10, 11, and 12. Both vortices 
of same sense. 

30°/0° 

3 

1.12 

Similar to flap/tip loading but with 
additional vortex of opposite sense 
shed from inboard edge of the flap. 
Simulates transport aircraft of 
aspect ratio 7 with inboard flap 
deflected by 30® and the outboard 
flap retracted. 

30°/30° 

5 

1.19 

Simulates transport aircraft of aspect 
ratio 7 with both flaps deflected by 30®. 

Elliptic 

8 

1.66 

Approximates, stepwise, an elliptic lift 
distribution. 

























Figure 4.- Experimental situation which allows for spanwise optical access. 
Typical of smaller facilities (e.g., Ames 7- by 10-Foot Wind Tmmel). 
Section lift coefficient can be determined from cc^ = u • di . 
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Figure 5.- Experimental situation which only allows for limited optical access 
Typical of large facility (Ames 40- by 80-Foot Wind Tunnel). Section lift 
coefficient determined by correcting measured cc^ using an equivalent 
vortex. 
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Z = Z2 


UNMEASURED 
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z = zi 


Figure 6.- Measured region determined by finite length of experimental traverse 
over only portion of wake survey line. Flow is modeled in unmeasured 
regions to generate a correction to measured value of 












Figure 7.- Flow chai-t for numerical wake survey experiment. 




2.0 





Figure 8.- Difference between "actual” and "measured" loading distribution for 
finite traverse distance. z = -z = 1.0, x = 0.9. 





Figure 9.- Flow chart for numerical correction to lift coefficient. 




Figure 10.- Corrected data showing coBiparison to ''actual" loading when lift 
coefficient = 1.08 is matched to within 1%. Equivalent vortex has 
strength of (cc„) = 1.42. 
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Figure 15.- Equivalent vortex correction applied to an approximately elliptic 

loading. 
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Figure 18.- Comparison of "actual" induced drag, -u (cc 

Z X, 

the distribution computed from equation (27). 30°/30 
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appendix a 


ANALYSIS OF THE SHEAR STRESS TERMS IN THE INTEGRAL MOMENTUM EQUATION 


To determine the relative importance of the shear forces on control 

surface (figs. 1 and 2) compared with the other terms 

assume that it is valid to conduct a two-dimensional analysis. The orders 
magnitude and the trends should he similar in three dimensions. 


MornGtitum Equation in Two Dimensions 

In the two-dimensional approximation of the flow, equation (3) simplifxes 



r 

^ J 

Lj^ thru 


(n • T)dil 


(Al) 


Separating into components, 


^ “ k 

pu u dz 
X z 

- / (pu ^ + p)dx 

- I pu u dz + 

A ^ ^ 

f (pu + p)dx 
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thr u Li^ 
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D = D 
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+"sf = j 

(pu ^ + p)dz - 
T ^ 

/ pu u dx - / 
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(pu^^ + 
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( ^ Vz 

dx + ^ (u 

• d£ 
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thru Lj^ 





Eliminating the static pressures with, p p,j, ^ 

nondimensionalizing , 
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^ ^ ^ A 3 ^ A^ 
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-2/ uu dx-f (u 
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2 - u ^)dz + 2 / u u dx 
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X z 
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(A3) 


Modeling of the Shear Stress Tensor 
The shear stress tensor is written in general form as 


9U; 


X 


T = y 


3x 


9u 9u 
X 4 . z 


9u 9u 
_x , z 

9z 9x 


9z 9x 
9u 

2 

9z 


If the flow is represented (fig. A-1) by a superposition of the free-stream 
flow, a bound vortex, and a viscous wake, 




- = (v + u” + A 

" X 2TT -^2 


/ w r x\- 

rj^ ^ \^y ■ Ztt 


where V is the strength of the bound vortex and the superscript w denotes 
a perturbation term due to the viscous wake. These wake terms vanish every- 
where except in the wake region. The stress tensor now becomes. 


9u 


w 


X 2r xz 


■e- 
T = yi 


9x TV 

„ w ' w 
9u 9u „ 2 

x . z , r x'" - z‘ 

+ — _ -I — 
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(A4) 


9z 


9x ir 


9z IT -r^ 


Since the trace of the matrix remains invariant and vanishes for an incompres- 
sible flow. 
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as expected. 
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9x 
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(A5) 


Evaluation of the Shear Forces 

To carry out the integration indicated by the last term in equations (A3)a 
we need to compute (x • ^ and (z • r). Using (A5), 


(x • t) = (1 0 )(?) = y 


3u 


X 2r xz 


9x 7T 


3u 3u ^ 

X i , z , r / x“^ - z‘ 
+ — :: — + 


3z 3x ir 


(z • ?) = (0 1 ) (t) = y 


30 ^^ 30 ^^ 

X . z , r 
+ — r 


3z 3x TT 

a w 

_ X . ^ xz 
3x TT ■r'+ 


(^)J 

i^j 






(A6) 




Equation (A 6 ) is valid everywhere. If, however, it is assumed that the 
shearing due to the viscous wake is only significant over a localized region 
along L^, then we have elsewhere. 


-2 


xz 


yr 

(x . t) = — 

" ' x2 ^ zZ 
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and (z • t) = 

TT 
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The shear terms in equations (A3) may now be written as 
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The integrations are greatly simplified and no generality is lost if we choose 
the symmetric controur in figure A-2. 


Using equations (A6) and (A7) in (A8), 



Since the integrands of the first and third terms are even in x and z, 
these terms cancel. Also, since the integrands of the second and last 
terms are odd in x and z, these terms cancel. Equation (A9) is analyzed 
similarly and the result is 
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thru 
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The last terms in (A3) may, therefore, be written nondimens ionally as 
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(5 . di = ^ 


thru Lj^ 


^ f 
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f /3u ^ 3u 

i, k * 

f 9u ^ 


dz 


thru 


where Re^ = pcV^/y is the flow Reynolds number based on wing chord. Clearly, 
when the velocity gradients in the wake are of order unity, then these terms 
are negligible for Re^ >> 1* This is certainly the case for most wind-tunnel 
flows . 
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WAKE DEFECTS 

WITH COMPONENTS 
w w 

U AND U 
X z 



Tigure A-1.- Flow modeled by superposition of free stream, bound vortex, 

and viscous wake. 


i—e, K ) 



(— e, -k) 


(e, — k) 


Figure A- 2.- Symmetrical integration contour. 




APPENDIX B 


EXPANSION OF THE SURFACE INTEGRALS IN EQUATIONS (9) 


The terms that are to be expanded are 


(u 


u )dx dz 


-2 / 

/g y z 

-2 I “ (u u )dx dz 

4 9y X y 


(Bl) 


Using the chain rule and then using the principle of continuity, (Bl) becomes 
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(B3) 


where the surface integrations are now shown explicitly by the double integral 

Evaluation of (B2) 

First term: integrate by parts to get 
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Third term: 
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Evaluation of (B3) 


(B4) 


First term: 
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Third term: 



(B3) now is 



Noting that 
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APPENDIX C 


THE INDUCED VELOCITY AT A POINT (x, y, z) DUE TO A HORSESHOE VORTEX 


To calculate the induced velocity at a point P, we use the Biot-Savart 
Induction Law 


-V 

j-5" r ds X r 

du = *7 — 

4tt „3 


(Cl) 


where the vector quantities are defined by ^gure C-1* The vector r from 
the source element of vorticity of length ds to point P is either 

'^+5 r_ , depending on the section of the horseshoe being considered. 

It can be shown that the induced velocity at P due to the bound vortex is 
given by 


• ^ X 

B Ait 


yn + y 
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Vo " y 


/x"' + (yQ + y)^ + /x^ + (y^ - y) 


+ z- 


^ (zx - xz) 


(C2) 


The induced velocity at P due to the trailing portion of the horseshoe on 
the positive y-side of the origin is given by 


4tt 


(y,, - y)’ 


+ 1 


+ z‘ 


[(yg - y)^ ■+ 2^1 ' [zy + (yg - y)2] 


(C3) 


and the induced velocity at P due to the trailing portion on the negative 
y-side of the origin is given by 


•V r r 

u = 
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+ 1 


^ + (yQ + y)^ + 2^ 


[^(Yq + y)^ + \[zy - (y-Q + y)z] (C4)^ 


Combining (C2) , (C3) , and (C4) and separating components. 
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A = 
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Nondimensionalizing u by V , and x, y, z by c, 


u = 1 + yAz 
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u = y(C - B)z 

y 


(C6) 


= -Y[Ax + B(yg - y) + C(y^ + y)] 

with A, B, and C unchanged but all terms now nondimensional, and y = cc^/Stt 
is the nondimensional vortex strength. 

To protect against singularities in the wake, a finite core diameter r^, 
is specified for the trailing portions of the horseshoe. This is accomplished 
by replacing 


B = 


r ^ L 
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X 


/x^ + (y^ - y) 
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+ 1 I if (y. “ y)2 + z2 < r 2 


C = 
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■^0 ^ c 


To allow for downward deflection of the wake, we assume that in a 
superposition of N horseshoes, each vortex trails downward at angle 
£j^(i = 1, 2 N) induced by the other vortices (see fig. C-2). 

We want to calculate the downwash velocity u„ induced on the bound 

i 

vortex at location y^ where Fj is shed. To do this we again use the 
Biot-Savart Induction Law. Assuming that all other vortices trail straight 


back, it can be shown that the induced velocity at yg due to the other 

i 

vortices on the (+) y-side of the centerline is given by 
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The induced velocity at y^ due to all vortices on the (-) y-side of the 
centerline is given by 

^ r 

1 

U = - — 

z 


y 

^ (>'0. ■*' yo.) 

3=1 \ 1 3/ 


(C8) 


Combining (C7) and (G8) , the total induced downwash is 
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If we take the induced angle as ^ ~V 
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•, then in nondimens ional form. 
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The flow field is now generated in the following manner: 

(i) Specify (x, y, z) where the velocity is desired. 

(ix) Determine rotated coordinates (x’ , y’ , z') for the ith horseshoe by 

cos e. 0 -sin e . 

X X 



0 


0 


sxn e. 0 cos e 

X X 



with computed by equation 
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(iii) 


! 


Compute velocities in the rotated coordinate system, (u ’ , u ' , 
’ ) , using equation (C6) . ^ 

(iv) Compute (u , u , u ) by rotating back to the x, y, z frame 
X y z 



Steps (i) through (iv) are carried out for each horseshoe and the total 
induced velocity at P is then obtained by linear superposition. 


48 




I 






I 




APPENDIX D 
DATA SUNtMARY 

Legend 

-f = "Measured" cc 
o = Corrected cc 
= "Actual" loading 

All data are for x = 0.9. 
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Figure D-5 
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